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Abstract: We study Kaluza-Klein (KK) graviton production in the large extra dimensions 
^H m °del via 2 jets plus missing transverse momentum signatures at the LHC. We make predictions 
^^for both the signal and the dominant Zjj and Wjj backgrounds, where we introduce missing 
QPy-dependent jet selection cuts that ensure the smallness of the 2-jet rate over the 1-jet rate. 

^With the same jet selection cuts, the distributions of the two jets and their correlation with 
*^>the missing transverse momentum provide additional evidence for the production of an invisible 

S_imassive object. 
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1. Introduction 

A general expectation in high energy physics today is that physics beyond the standard model 
(BSM) should emerge at TeV energies. This belief is founded on the observation that the elec- 
troweak symmetry breaking scale of the SM cannot be made stable against quantum corrections 
without invoking new physics at the TeV scale. With this in mind, an enormous international 
effort is being poured into the construction of the 14- TeV Large Hadron Collider (LHC) at 
CERN. Apart from supersymmetry (SUSY), models with large extra space dimensions, such as 
the one proposed by Arkani-Hamed, Dimopoulos, and Dvali (ADD) [Q], provide an alternative 
possibility in this direction that is most exciting. 

In the D = 4+5 dimensional ADD model, the SM particles live in the usual 3+1— dimensional 
space, while gravity can propagate into the additional ^-dimensional space, which is assumed 
for simplicity to be compactified on the <5-dimensional torus T s with a common radius R. Then 
the 4-dimensional Planck scale Mpi is related to the fundamental scale M s as follows 

M Pl = 8irR s M!> +2 , (1.1) 



where M s ~ TeV is possible for large compactification radius R. According to Eq. (LI), one 
can expect that deviations from the usual Newtonian gravitational force law will appear at a 
distance around R ~ 0.83 x 10~ 16+ ~mm(2.4TeV/M s ) 1+ « . Terrestrial experiments gave the 
limit R < 0.2 mm by probing gravitational forces directly 0. For 5 = 2 this translates into 
Ms ^1-5 TeV, while for 5 > 2, there are no strong limits on M s . 

In our four dimensional space-time, there appear Kaluza-Klein (KK) towers of massive spin- 
2 gravitons in the ADD model which interact with the SM fields. The interaction Lagrangian is 
given by §, § 

£ int = - * ^G^T^, (1.2) 
M Pl H 

where the massive gravitons are labeled by a 5-dimensional vector of positive integers, ri = 
(ni,n2, ■•,",5), Mpi = Mpij\f%Tx ~ 2.4 x 10 18 GeV is the reduced four dimensional Planck scale, 
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and T au is the energy-momentum tensor of the scattering fields. The n-th KK mode graviton 
mass squared is = \n\ 2 /R 2 . From Eq. ( |1.2| ) one can derive the relevant Feynman rules, 
some of which can be found in Ref. [||, ||. For M s = 1 TeV and 5 =4, 6 and 8, the mass gap of 
the KK modes is Am = R^ 1 ~ 20 keV, 7 MeV and 0.1 GeV 0|, respectively. Thus the spectrum 
of KK modes can be treated as continuous for 5 < 6, and approximated by the mass density 
function Q 

p(m) = Ss-ij^-m 5 - 1 , with S S -i = T ( S / 2 ) • ( L3 ) 

So far the strongest constraints for 5 < 4 extra dimensions come from astrophysics and 
cosmology, but they can be relaxed and do not diminish the importance of collider phenomenol- 
ogy ||. We therefore discuss the 5 = 3 case as well. In collider experiments, there are two 
classes of effects that can probe large extra dimensions: virtual KK tower exchange between the 
SM particles and a real graviton emission. Since the couplings of the graviton with matter are 
suppressed by inverse power of Mpi, graviton direct production leads to missing energy signals. 
Detailed studies have been performed at LEP || and the Tevatron J7|, by searching for gravitons 
in the processes e + e~ — » j(Z) + E miss and pp —* jQet) + P™ 1SS . The combined LEP 95% CL 
limits are M s > 1.60, 1.20, 0.94, 0.77, 0.66 TeV for 5 = 2,- • • ,6 respectively, and the Tevatron 
95% CL limits are M s > 1.18, 0.99, 0.91, 0.86, 0.83 TeV. For the LHC, graviton production with 
a monojet has been investigated in detail and found to have strong ability to probe up to much 
higher extra dimension scale Q. There is very little information on the underlying physics, 
however, in the transverse momentum and the rapidity of the single jet. One may thus wonder 
whether additional jets in graviton production can be used as a more sophisticated probe. 

In this paper, we study graviton production with 2 jets at the LHC in large extra dimen- 
sions, in comparison with the dominant Zjj background ||], and examine if the 2-jet rate and 
correlations can give us more information about the mass scale of the missing object, in addition 
to the missing Pj- distribution. At the same time we address the questions to what extent high 
transverse momentum graviton production will indeed emerge as a monojet signature. We will 
show that higher order QCD effects do lead to a more complex event structure, with multiple 
jets in the 100 GeV Pt range. The remainder of this work is organized as follows: In Section II 
we present the calculations. In Section III we give numerical results and discussions. Section IV 
contains our conclusions. 



2. Calculations 

We are considering the QCD production of a graviton with 2 jets at the LHC, pp — > jjG n , 
including all the possible subprocesses, among which gg — > ggG n , gq — > gqG n and qq( •* — ► qq( ^G n 
play the most important role. Representative Feynman diagrams are shown in Fig. |]. In addition 
to the QCD processes of Fig. 1, we have also calculated the electroweak (EW) contributions to 
jjG n production. In particular, we have determined the graviton production cross sections from 
weak boson fusion (WBF) processes. However, the WBF cross sections for jjG n production 
represents a small correction, which is below 1%, even when imposing typical cuts to enhance 
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WBF over QCD sources [ 10 1 . Thus, WBF processes do not appear as a promising avenue for 
studying graviton production at the LHC, and we do not include them in the results below. 

Significant background can come from any processes leading to two jets and missing trans- 
verse momentum, among which we consider the most important one, namely Zjj production 
with subsequent decay Z — > vv |J. We also studied another class of processes which could be 
significant at least when missing Pj- is not too large. This can arise from QCD production of 
Wjj with subsequent decay — > l^v when the charged leptons I = e, /x, r are not identified. 
Here we follow the procedure of Ref. [Til] . 




q 

Gn 



Gn / 




Figure 1: Representative Feynman diagrams for the subprocesses (a) qq^ — > qq^G n , (b) qg — ► qgG„ 



and (c) gg — > ggG n which contribute to the dijet plus graviton production process, pp — > jjG r , 
The gravitons arc emitted from each of the solid points in the diagrams. 
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The signal and background are simulated at the parton level with full tree level matrix 
elements. The amplitudes are calculated by the helicity amplitude technique and we have 
added all the relevant HELAS subroutines |T3] for the massive graviton and its interactions based 



on the effective Lagrangian of Eq. (1.2). For the background, we used the codes based on Ref. jR 
and checked by MadGraph/MadEvent |14|j . For the signal, we have performed two independent 
calculations to check each other. Firstly, we wrote our simulation codes based on the ones 
generated for calculations of graviton radiation at linear colliders jl5), which matches closely 
with similar other calculations |fj|. Secondly, we have also implemented ADD spin-2 gravitons 
into MadGraph/MadEvent. We find agreement between the two independent calculations. 

Additional checks were carried out for the Born level amplitude. The most useful check is 
provided by the Ward identities arising from general coordinate invariance, which constitutes an 
essential feature of any theory involving gravity. We can write the amplitude for the emission 
of any graviton in the form 

A n (k, Pi ) = T^(k, Pi )e<$*(k) , (2.1) 



where pi are the momenta of the external SM particles, e^J(k) is the polarization tensor for 
the n-th (massive) graviton mode with its momentum k. The tensor T^ y {k,pi) is the same 
for all the graviton modes, including the massless mode e^}{k), which is the graviton of the 
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four-dimensional Einstein gravity. This must now satisfy the Ward identities 

k^T^pi) = hTT^fapi) = , (2.2) 

where we note that 

N 

T^(k, Pl )=Y,T^(k,p t ) (2.3) 
i=i 

with j indicating the j-th. diagram, as above. The consistency check therefore requires a perfect 
cancellation between Feynman graphs, for each choice of p or v, which is highly sensitive to 
errors in signs and factors. Our numerical check in Ward identities confirms the cancellation up 
to the expected accuracy of our numerical programs. 

Once dealing with this effective low-energy theory, one concern is its behavior above the 
ADD fundamental scale (M s ). One unitarity criterion which we have implemented is that the 
tower of gravitons being produced does not extend in mass beyond the ADD scale (Mg„ < M 8 ). 
Moreover, we will also present the results with a hard truncation scheme, by setting the cut 
^truncation < M s , where the truncation parameter is set as the root of the partonic center-of- 
mass energy, 

^truncation — V^§, (^4) 

which is a quite conservative choice. 
3. Results and Discussions 

In the tree level numerical calculations, we identify massless partons with jets which must satisfy 
the angular cuts 

ARjj = VAr? 2 + A0 2 > 0.7 , \rjj\ < 4.5. (3.1) 

Here rj is the pseudorapidity of the jets and <f> is the azimuthal angle around the beam direction. 
Unless specified otherwise we further require 

> 6 GeV x Jpfssjl GeV, (3.2) 

P^ iss > 1 TeV . (3.3) 

We employ CTEQ6L1 parton distribution functions (PDF) 0] throughout, with the factoriza- 
tion scale chosen as \if = min(Pr) of the jets which satisfy the above cuts. The QCD coupling 

is set to the geometric mean value, a s = \J a s (P^ 1 ) a s (Pj, 2 ) . For the ADD parameter, we first 
focus on the 5 = 4 and M s = 5 TeV case in Figs. ^-||, and then discuss the ADD scale sensitivity 
and present the differential distributions for 5 = 3, 4, 5, 6 cases, respectively. For 5 > 6 large 
extra dimensions, the total cross sections for real graviton production are smaller and not easy 
to detect at the LHC, thus we will not discuss them here. 

Notice that throughout our calculations, we follow the notation of Ref. || which differs 
from the one in Ref. 0] mainly by a different factor in the relation between R and M s in (4+5)- 
dimensional space. Though this factor is crucial in comparing results and quantifying discovery 
potentials, one can simply convert results from one notation to the other by multiplying a 5 - 
dependent factor. 
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Figure 2: cut dependence of the dijet cross sections for signal (a) and background (b) at the LHC 
in various missing transverse momentum bins when applying the cuts of Eq. (3.1). The open circles show 
the monojet cross section in the same missing Pp bin. 
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The P T cut of Eq. Q3,2| ) is chosen to in- 
sure perturbative ordering of the tree level 
cross sections throughout phase space, i.e. 
we want to keep the dijet cross sections be- 
low the corresponding 1-jet inclusive cross 
sections. In order to motivate our choice, we 
show, in Fig. |2|, the P T cut dependence of the 
total cross sections for the dijet plus missing 
Pt events at the LHC. Fig. 2(a) is for the 
signal process pp — > jjG n X, and Fig. 2(b) 
is for the background process pp — » jj(Z — > 
vv}X. Each line shows the dijet cross sec- 
tion for the missing Pt in a 100 GeV bin 
between 200 GeV and 1400 GeV. The open 
circle along the lines shows the correspond- 
ing monojet cross section in the same miss- 
ing Pt bin. In smaller bins, the jjG n 
total cross sections drop faster, with increas- 
ing P T cut, than the background ones, which is due to the soft and collinear Z boson emission 
along one of the QCD jet directions in the Zjj background. 

The P™ m values at open circle points in Figs. 2(a) and (b) tell the jet Pt threshold below 
which the dijet cross section is larger than the monojet cross section, and hence our perturbative 
results cannot be trusted. Below these threshold values one should expect multiple soft jet 
emission to appear. For a missing Pt of 1 TeV, for example, gluons with Pt ^$ 140 GeV are in 
the soft range, and several such "soft" gluon jets are expected in a typical graviton or background 



1000 

(GeV) 



Figure 3: Missing transverse momentum depen- 
dence of the P™* value of equal 2-jet and 1-jet cross 



sections. Our jet selection cut Eq. (3^) is also pre- 
sented. 
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event. Since these gluons are readily observable as distinct jets in the experiment, an actual 
monojet event with missing transverse momentum in the TeV range and no additional jets with 
Pt ~ 30 GeV, is a very rare event. 

To get results which are perturbatively 
reliable, we introduce the missing Pt depen- 
dent jet selection cut of Eq. (|3.2|) such that 
the dijet to monojet cross section ratio is 
always smaller than unity, while keeping as 
many dijet events as possible. We show in 
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Fig. 3 our jet Pt cut of Eq. (3^), together 
with the Pt threshold values for the sig- 
nal (dashed red) and the background (dot- 
ted blue) above which the dijet cross sec- 
tion is smaller than the monojet one. We 
find from this figure that a larger fraction 
of graviton dijet events survives the cut for 
P^ iss < 900 GeV, while Z + jet(s) events 
obtain higher jet multiplicities at higher P™ 1SS . 
This is partly because of the higher hard 
scattering scale of the graviton events at small 

Pj? lss , and partly because of the importance of collinear Z boson emission in high P™ 1SS back- 
ground events. 
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Figure 4: P™ lss dependence of the total cross sec- 
tions for the signal and background when applying the 
cuts of Eqs. 
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Figure 5: The dependence on ADD scale M s of the total cross sections for jG n (a) and jjG n (b) 
productions at the LHC. The SM background results also have been plotted. The dashed lines are for 



Ma„ < M s , while the solid lines are for y/I < M s . Event selection criteria chosen as in Eqs. (3.1) - (f3~3|) 



In Fig. ||, we show the P™ 1SS spectrum for both the 1-jet and 2-jets signal and background 
processes. One finds that the > 1 TeV requirement in Eq. ( |3.3| ) reduces the background 
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sufficiently. Also included is the contribution coming from the Wjj background which falls 
sharply for large P™ 1SS and which yields an additional background contribution of 0.15 fb above 
P™ 1SS = 1 TeV, i.e. it is negligible compared the the Zjj background. A further improvement 
of the signal to background ratio is possible by tightening the P™ 1SS cut, but this will not be 
pursued in the following. 




Figure 6: Distributions of the azimuthal angle separation between the two jets (left), and the min- 
imal azimuthal angle separation between the two jets and missing transverse momentum (right), for 
background and graviton signal. 

In Fig. [5|, we present the ADD scale de- 
pendence of the total cross sections for jG n 
and jjG n production, for the 5 = 3,4,5,6 
cases. The SM juu and jjvv background 
results are also plotted. Our results for jG n 
production with the hard truncation scheme 
agree with the results in Ref. || within about 
5 percent, which may be due to different 
PDF and scale choices. However, the re- 
sults without truncation have larger differ- 
ences especially at small M s , because the 
unitarity criterion Ma„ < M s is used in 
our paper as mentioned above, while not in 
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Table 1: Maximum ADD scale M s sensitivity which 
can be reached by studying the 2-jet (1-jet) and miss- 
ing transverse momentum signal at the LHC, with in- 
tegrated luminosity C =100 fb -1 , assuming the sys- 
tematic error to be 10%. 



Ref. |Q. We have also performed the same sensitivity analysis as in Ref. ||, considering the 
integrated luminosity C = 100 fb -1 , where the systematic error in the background (assumed to 
be 10%) dominates over the statistical error. The sensitivity range is defined by 

VjjGMjGn) > 5 x 10% x o- background = 1.93 (2.45) fb. (3.4) 

The resulting max M s sensitivity results are shown in table. I. The 2-jet sensitivity is only 
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slightly lower than for the 1-jet case. Moreover, the larger 5 is, the sooner the non-perturbative 
region is reached, thus the larger is the difference between max M s sensitivities in no truncation 
and hard truncation cases. 

The most distinct difference between the signal and background in the dijet plus missing 
Pt events is found in azimuthal angle correlations between the 2 jets and the missing trans- 



verse momentum. The 4>jj and IXLin{ (D ■ pmiss ) distributions with the cuts (|3.1| )-( |3~3| ) are shown 
in Figs. ||, for M s = 5 TeV and 5 = 3,4, 5, 6, respectively. The Zjj background shows a clear 
enhancement for back to back jets, reflecting collinear Z emission along the direction of one of 
the jets. Due to the heavier masses of the typical graviton KK modes, such collinear "jet frag- 
mentation" contributions are absent for the signal. This significant difference of the azimuthal 
angle distributions can provide a powerful tool to test for heavy graviton emission: the relative 
suppression of the jet fragmentation contribution in the data would be a direct sign for a very 
massive object as a source of the missing transverse momentum. 

Finally, we would like to comment on 
the ratio of the 2-jet over the 1-jet rate. 
Since the typical graviton KK modes are 
much heavier than the Z boson, thus pro- 
viding for a much harder event, one would 
expect naively that the ratio for the signal 
should always be larger than the one for the 
background. However, that is not the case 
as can be seen from the results in Fig. ^, es- 
pecially at large P™ 1SS . The reason can be 
traced back to final state parton emission, 
which is more prominent for the SM back- 
ground. Based on this observation, we show, 
in Fig. [?], the 2-jet over 1-jet ratio for sig- 
nal and background clS cl function of Pj? lss . 
The upper two curves are from the results 
in Fig. 0. For the lower two curves, we set 
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Figure 7: The 2-jet over 1-jet ratio for signal and 
background as a function of P™ ISS , with or without the 
cut |A?7jj| > 2 and \<j>jj — 7r| > 0.7. Additional cuts of 



Eqs. (3.1), (3.2) applied 



additional cuts lA^-j > 2 and \4>jj — ir\ > 

0.7, in order to reduce final state parton emissions. One now finds a higher dijet fraction for the 
graviton events, largely from initial state radiation, reflecting the harder collision scale of the 
graviton emission events at all missing Py. 



4. Summary 

The production of stable Kaluza-Klein gravitons at the LHC can be observed when they are 
produced at large transverse momentum, giving rise to a large missing Pp signature. For 
of order 1 TeV or larger, the signature will rarely be a monojet signal, however, because multiple 
"soft" gluon emission will produce events with several jets balancing the transverse momentum 
of the graviton. We have calculated the order a 2 s graviton plus dijet, jjG n , cross section at the 
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LHC and find that it saturates the leading order monojet cross section for additional "soft" jet 
Pt in the 100 to 150 GeV range, thus establishing the typical scale for multiple jet emission. 

The Pfi 1SS distribution of the jjG n cross section is strongly influenced by the minimal jet 
Pt requirements imposed on the partons. We note that defining the dijet cross section with a 
typical constant jet Pt cut, independent of the hardness of the event, will invariably lead to 
the cross section not being trustworthy at sufficiently high P™ 1SS , where saturation of the LO 
monojet cross section happens at higher jet Pt already. Not considering this effect, via a sliding 
jet Pt cut as in Eq. (|3.2| ), may produce an unrealistically hard missing transverse momentum 
distribution, which can lead to an overestimate of the LHC sensitivity to graviton production. 

In addition to the order a 2 s process, jjG n production via gluon exchange, we have also 
calculated the corresponding electroweak process qq — > qqG n via weak boson fusion. However, 
the cross section for the latter is always strongly suppressed. Even with typical weak boson 
fusion cuts we have been unable to find phase space regions where the electroweak process 
contributes to overall jjG n production at more than the percent level, while also being visible 
above the SM background. We conclude that weak boson fusion is not a promising process for 
Kaluza-Klein graviton production at the LHC. 

The multijet characteristics described above for the signal are also expected for the dominant 
SM background, Zjj production with subsequent decay of the Z boson to neutrino pairs. The 
multijet features are simply a reflection of the hardness of the event, as specified by the large 
missing transverse momentum. We have found one feature, however, which distinguishes signal 
and background processes. For missing Pt in the TeV range the Z mass becomes negligible 
and jet fragmentation into a collinear Z becomes an important part of the SM background. 
This contribution is most readily seen in the azimuthal angle correlations of the jets, with 
a sizable fraction of nearly back-to-back dijet events. The large average mass of the produced 
gravitons strongly suppresses such a contribution for the signal. The resulting distinctly different 
azimuthal angle distributions may be useful to verify that an observed excess of high P™ 1SS events 
is indeed due to the production of an invisible very massive particle. 
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